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[i]  A  new  concept,  valid  prediction  period  (VPP),  is  presented  here  to  evaluate  ocean 
(or  atmospheric)  model  predictability.  VPP  is  defined  as  the  time  period  when  the  prediction  error 
first  exceeds  a  predetermined  criterion  (i.e.,  the  tolerance  level).  It  depends  not  only  on  the 
instantaneous  error  growth  but  also  on  the  noise  level,  the  initial  error,  and  the  tolerance  level.  The 
model  predictability  skill  is  then  represented  by  a  single  scalar,  VPP.  The  longer  the  VPP,  the  higher 
the  model  predictability  skill  is.  A  theoretical  framework  on  the  basis  of  the  backward  Fokker- 
Planck  equation  is  developed  to  determine  the  mean  and  variance  of  VPP.  A  one-dimensional 
stochastic  dynamical  system  [Nicolis,  1992]  is  taken  as  an  example  to  illustrate  the  benefits  of  using 
VPP  for  model  evaluation.  Index  Terms:  4263  Oceanography:  General:  Ocean  prediction; 
4255  Oceanography:  General:  Numerical  modeling;  3367  Meteorology  and  Atmospheric 
Dynamics:  Theoretical  modeling;  KEYWORDS:  backward  Fokker-Planck  equation,  instantaneous 
error,  Lorenz  system,  predictability,  tolerance  level,  valid  prediction  period 


1.  Introduction 

[2]  A  practical  question  is  commonly  asked:  How  long  is  an 
ocean  (or  atmospheric)  model  valid  once  being  integrated  from  its 
initial  state,  or  what  is  the  model  valid  prediction  period  (VPP)?  To 
answer  this  question,  uncertainty  in  ocean  (or  atmospheric)  pre¬ 
diction  should  be  investigated.  It  is  widely  recognized  that  the 
uncertainty  can  be  traced  back  to  three  factors  [Lorenz,  1969,  1984]: 

( 1 )  measurement  errors,  (2)  model  errors  such  as  discretization  and 
uncertain  model  parameters,  and  (3)  chaotic  dynamics.  The  meas¬ 
urement  errors  cause  uncertainty  in  initial  and/or  boundary  con¬ 
ditions  [e.g Jiang  and Malanotte-Rizzoli,  1999],  The  discretization 
causes  small-scale  “subgrid”  processes  to  be  either  discarded  or 
parameterized.  The  chaotic  dynamics  may  trigger  a  subsequent 
amplification  of  small  errors  through  a  complex  response. 

[3]  The  three  factors  cause  prediction  error.  For  example,  an 
experiment  on  the  Lorenz  system  was  recently  performed  [Chu, 
1999]  through  perturbing  initial  and  boundary  conditions  by  the 
same  small  relative  error  (1(L4).  The  vertical  boundary  condition 
error  is  transferred  into  the  parameter  error  after  turning  the 
Saltzman  model  [Saltzman,  1962]  into  the  Lorenz  model  [Lorenz, 
1963]  model  using  a  variable  transform.  The  relative  error  is 
defined  by  the  ratio  between  the  rms  error  and  rms  of  the  three 
components.  The  Lorenz  system  has  a  growing  period  and  an 
oscillation  period.  With  the  standard  parameter  values  as  used  by 
Lorenz  [1963],  the  growing  period  takes  place  as  the  nondimen- 
sional  time  from  0  to  22;  and  the  oscillation  period  occurs  as  the 
nondimensional  time  from  22.  During  the  growing  period  the 
relative  error  increases  from  0  to  an  evident  value  larger  than  1 
for  both  initial  and  boundary  uncertainties.  During  the  oscillation 
period  the  relative  error  oscillates  between  two  evident  values:  4.5 
and  0.1  for  the  initial  uncertainty  and  5.0  and  0.2  for  the  boundary 
uncertainty. 

[4]  Currently,  some  timescale  (e.g.,  e-folding  scale)  is  computed 
from  the  instantaneous  error  (defined  as  the  difference  between  the 
prediction  and  reality)  growth  to  represent  the  model  VPR  Using 
instantaneous  error  (IE),  model  evaluation  becomes  stability  anal- 
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ysis  on  small-amplitude  errors  in  terms  of  either  the  leading 
(largest)  Lyapunov  exponent  [e.g.,  Lorenz,  1969]  or  calculated 
from  the  leading  singular  vectors  [e.g.,  Farrell  and  Ioannou, 
1996a,  1996b],  The  faster  the  IE  grows,  the  shorter  the  e-folding 
scale  is  and,  in  turn,  the  shorter  the  VPP  is. 

[5]  For  finite  amplitude  IE,  however,  the  linear  stability  analysis 
becomes  invalid.  The  statistical  analysis  of  the  IE  (both  small- 
amplitude  and  finite  amplitude)  growth  [Ehrendorfer,  1994a, 
1994b;  Nicolis,  1992],  the  information  theoretical  principles  for 
the  predictability  power  [Schneider  and  Griffies,  1999],  and  the 
ensembles  for  forecast  skill  identification  [Toth  et  al.,  2001] 
become  useful. 

[6]  The  probabilistic  properties  of  IE  are  described  using  the 
probability  density  function  (PDF)  satisfying  the  Liouville  equa¬ 
tion  or  the  Fokker-Plank  equation.  Nicolis  [1992]  investigated  the 
properties  of  the  IE  growth  using  a  simple  low-order  model 
(projection  of  Lorenz  system  into  most  unstable  manifold)  with 
stochastic  forcing.  A  large  number  of  numerical  experiments  were 
performed  to  assess  the  relative  importance  of  average  and  random 
elements  in  the  IE  growth. 

[7]  Although  familiar  and  well  understood,  the  IE  growth  rate 
is  not  the  only  factor  to  determine  VPP.  Other  factors,  such  as  the 
initial  error  and  tolerance  level  of  prediction  error,  should  also  be 
considered.  The  tolerance  level  of  prediction  error  is  defined  as  the 
maximum  allowable  forecast  error.  For  the  same  IE  growth  rate 
the  higher  the  tolerance  level  or  the  smaller  the  initial  error,  the 
longer  the  model  VPP  is.  The  lower  the  tolerance  level  or  the 
larger  the  initial  error,  the  shorter  the  model  VPP  is.  Thus 
the  model  VPP  is  defined  as  the  time  period  when  the  prediction 
error  first  exceeds  a  predetermined  criterion  (i.e.,  the  tolerance 
level  e).  The  model  predictability  is  then  represented  by  a  scalar 
VPP.  The  longer  the  VPP,  the  higher  the  model  predictability  is.  In 
this  study  we  develop  a  theoretical  framework  for  model  predict¬ 
ability  evaluation  using  VPP  and  illustrate  the  usefulness  and 
special  features  of  VPP.  The  outline  of  this  paper  is  depicted  as 
follows.  Description  of  prediction  error  of  deterministic  and 
stochastic  models  is  given  in  section  2.  Estimate  of  VPP  is  given 
in  section  3.  Determination  of  VPP  for  a  one-dimensional  sto¬ 
chastic  dynamic  system  is  discussed  in  section  4.  The  conclusions 
are  presented  in  section  5. 
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2.  Prediction  Error 

2.1.  Dynamic  Law 

[8]  Let  x(l)  =  [x(1)(0,  x<2)(f),  . . .,  x(n>(t)\  be  the  full  set  of 
variables  characterizing  the  dynamics  of  the  ocean  (or  atmosphere) 
in  a  certain  level  of  description.  Let  the  dynamic  law  be  given  by 

f  =  f(M),  (1) 

where  F  is  a  functional.  Detenninistic  (oceanic  or  atmospheric) 
prediction  is  to  find  the  solution  of  equation  (1)  with  an  initial 
condition 


with  the  IE  growth  rate  and  the  corresponding  e-folding  timescale 
as  the  measures  of  the  model  predictability  skill.  It  is  reasonable  to 
assume  linear  error  dynamics  for  small-amplitude  errors, 


f(z,  t )  =  A(t)z. 


The  first  Lyapunov  exponent  is  defined  as 


X  =  lim  sup 

/— >  OO 


ln(||$(t,to)||) 


(9) 


(10) 


where  4>(t,  t0)  is  calculated  by  \Coddington  and  Levinson ,  1955] 


x(t0)  =  Xo, 


(2) 


where  x0  is  an  initial  value  of  x. 

[9]  With  a  linear  stochastic  forcing  q{t)x,  equation  (1)  becomes 

dx 

-Jt  =  f(x,  t)  +  q(t)  (x).  (3) 


<£(?,  fo)=I  +  J  A (s)ds  +  j  A (r)dr  J  A (s)ds  +  . . . , 

to  to  to 


where  I  is  the  unit  matrix.  Usually,  the  e-folding  scale  relating  to 
the  IE  growth  rate  (or  the  Lyapunov  exponent)  is  used  to  represent 
VPP. 


Here  q(t)  is  assumed  to  be  a  random  variable  with  zero  mean, 

<?(0>  =  °>  (4) 

and  pulse-type  variance 

{q(t)q{t'))  =q2b(t-t'),  (5) 

where  the  bracket  <  >  is  defined  as  the  ensemble  mean  over 
realizations  generated  by  the  stochastic  forcing,  5  is  the  Delta 
function,  and  q2  is  the  intensity  of  the  stochastic  forcing. 

2.2.  Model  Error 

[10]  Let  y(f)  =  [y(1)(f),  y<2)(0>  ■  ■  •,  y("\t)\  be  the  estimate  of  x(t) 
using  the  prediction  model  equation  (1)  or  equation  (3)  with  an 
initial  condition 


3.2.  Direct  Calculation 

[12]  When  state  errors  grow  to  finite  amplitudes,  the  linear 
assumption  equation  (9)  is  no  longer  applicable,  and  the  nonlinear 
effect  should  be  considered.  However,  the  prediction  error  cannot 
be  less  than  the  noise  level  (minimum  limit)  f(noise  and  greater  than 
the  tolerance  level  (maximum  limit)  e.  The  ratio  between  the 
maximum  and  minimum  limits  is  usually  large: 

D  =  e/tnoi*  »  !»  (n) 

in  ocean  models.  Thus  the  rms  error  of  prediction  is  bounded  by 
the  two  limits 

incise  <  A*)  <  E2-  (12) 


y(*o)  =  y0- 


(6)  Two  ellipsoids,  SE  and,  5^  are  defined  by 


The  prediction  error  vector  is  defined  by  T.  .  ->  r/  , 

•2(z)  =  E-  -Z(z)  =  ^,noise  , 


Z  =  X  y, 

at  any  time  t  (>t0),  and  the  initial  error  vector  is  defined  by 


Zo  =  Xo  -  y0. 


If  the  components  [x(1,(f),  x(2\t),  . . x<n>(t)\  are  not  equally 
important  in  terms  of  prediction,  the  uncertainty  of  model 
prediction  can  be  measured  by  the  rms  error 

A  z)  =  (z'Wz),  (7) 

where  W  is  the  diagonal  weight  matrix,  the  superscript  t  denotes 
the  transpose  operator,  andthe  bracket  represents  the  ensemble 
average  over  realizations  generated  by  stochastic  forcing,  uncertain 
initial  conditions,  and  uncertain  model  parameters. 


with  y (t)  as  the  center  (Figure  1). 

[13]  VPP,  represented  by  a  time  period  ( t  —  t0)  at  which  z  (the 
error)  goes  out  of  the  ellipsoid  Se(t),  is  a  random  variable,  whose 
conditional  PDF  P[(f  —  fo)|zo]  satisfies  the  backward  Fokker- 
Planck  equation  [Pontryagin  et  al.  1962\  Gardiner,  1983;  Ivanov 
et  al. ,  1994] 


dP 

~dt 


[f(zo,0] 


dP 

dz0 


1  2  2  92P  n 

2  C/ZqOZq 


(13) 


To  solve  equation  (13),  one  initial  condition  and  two  boundary 
conditions  (with  respect  to  z0)  are  needed.  Since  the  initial  error 
z0  is  always  less  than  the  given  tolerance  level  (i.e.,  always 
inside  the  ellipsoid  SE(i)),  the  conditional  PDF  of  VPP  at  t0  is 
given  by 


3.  VPP 

3.1.  Indirect  Estimate  From  IE  Growth  Rate 

[ll]  Traditionally,  the  stability  analysis  is  used  to  investigate  the 
small-amplitude  error  dynamics  due  to  the  initial  condition  error. 


P[(0)|z0]  =  1.  (14) 

If  the  initial  error  z0  reache  the  tolerance  level  (i.e.,  z0  hits  the 
boundary  of  SE(t0)),  the  model  loses  prediction  capability  initially: 


z(to)  =  z0, 


t  =  to, 


(8) 


P[(t  -  t0)  |  z0]  =  0  at  J (z0)  =  £2, 


(15a) 
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Figure  1.  Phase  space  trajectories  of  model  prediction  y  (solid 
curve)  and  reality  x  (dotted  curve)  and  error  ellipsoid  Se(t)  centered 
at  y.  The  positions  of  reality  and  prediction  trajectories  at  time 
instances  are  denoted  by  asterisks  and  open  circles,  respectively.  A 
valid  prediction  is  represented  by  a  time  period  ( t  —  t0)  at  which  the 


the  PDF  of  VPP  still  varies  with  time  (following  the  backward 
Fokker-Planck  equation,  equation  (13)): 


dP 

Ht 


f(zo) 


dP 

dz0 


1  2  2  d2P 

2  ^  0  dz0dz0 


=  0. 


We  multiply  this  equation  by  ( t  —  t0)  and  (f  —  t0)2,  then  integrate 
with  respect  to  t  from  t0  to  oo  and  obtain  the  mean  VPP  equation 


rt  ,  ?2Zo 

f(z0)T2— +  »  =~1 

C/Zq  2  OZqOZq 


(18) 


and  the  VPP  variability  equation 


<fz |  d2T2 

2  dz0dzo 


— 2ti  . 


(19) 


Here  the  expression 


OO 

J  P[(t  -  to)  I  Zo \dt  =  1 

to 


is  used.  Both  equations  (18)  and  (19)  are  linear,  time-independent, 
and  second-order  differential  equations  with  the  initial  error  z0  as 
the  only  independent  variable.  Two  boundary  conditions  for  iq  and 
t2  can  be  derived  from  equations  (15a)  and  (15b), 


which  is  the  absorbing-type  boundary  condition.  If  the  initial  error 
reaches  the  noise  level  (i.e.,  z0  hits  the  boundary  of  Sip(to)),  the 
boundary  condition  becomes  [ Gardiner ,  1983] 


dP[(t  -  t0)  |  zo] 


=  0,  J(z0)  =  Q 


2 

■noise  ’ 


(15b) 


which  is  the  reflecting  boundary  condition. 

[14]  The  Ath  moment  (A  =  1,  2,  . . .)  of  VPP  is  calculated  using 
PDF: 


Tl  =  0, 

T2  =  0, 

J{  zo)  =  e2; 

(20) 

<9ti 

i9t2 

— -  =  0 
dz0  ’ 

— -  =  0 
dzo  ’ 

OJ 

’S 

H 

o 

N 

(21) 

Example 

4.1.  One-Dimensional  Stochastic  Dynamical  System 

[16]  We  use  a  one-dimensional  probabilistic  error  growth  model 
[Nicolis,  1992] 


T k  (zo)  =  k  I  P[(t  -  to)  I  z0]  (f  -t0)k  1  dt,  k  =  1, . . . ,  00. 

to 

(16) 

The  mean  and  variance  of  VPP  can  be  calculated  from  the  first  two 
moments 


-jt  =  (oi  -  gi2)  +  v(t)i,  0  <  £,  <  oo,  (22) 

as  an  example  to  illustrate  the  procedure  in  computing  mean  VPP 
and  VPP  variability.  Here  the  variable  £,  corresponds  to  the  positive 
Lyapunov  exponent  ct,  g  is  the  nonnegative  parameter  whose 
properties  depend  on  the  underlying  attractor,  and  v(t)(,  is  the 
stochastic  forcing  satisfying  the  condition 


(t)  =  Tl 


(17a)  (v(t))  =  0,  (v(t)v(t'))  =  q2i{t  ~  !)■ 


(6t2)  =  t2  -  t2, 


(17b) 


where  the  bracket  denotes  the  average  overrealizations. 


3.3.  Autonomous  Dynamical  System 

[is]  The  predictability  is  usually  time-dependent  in  ocean  (or 
atmospheric)  systems  [Toth  el  al,  2001],  Even  for  an  autonomous 
dynamical  system, 


Without  the  stochastic  forcing  v(f)(,  the  model  (22)  becomes  the 
projection  of  the  Lorenz  attractor  onto  the  unstable  manifold. 

4.2.  Equations  for  the  Mean  and  Variance  of  VPP 

[17]  How  long  is  the  model  (22)  valid  once  being  integrated 
from  the  initial  state?  Or  what  are  the  mean  and  variance  of  VPP  of 
equation  (22)?  To  answer  these  questions,  we  should  first  find  the 
equations  depicting  the  mean  and  variance  of  VPP  for  equation 
(22).  Applying  the  theory  described  in  sections  3.2  and  3.3  to  the 
model  (22),  the  backward  Fokker-Planck  equation  becomes 


dP 

Hi 


IHo  -  g£o] 


dP  1  2  d2P 
dto~2q  Wo 


=  0, 


f  =  f(z„), 


(23) 


X  -  4 
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Figure  2.  Contour  plots  ofTi(£0,  £noise,  e)  versus  (£0,  £noise)  for  four  different  values  of  e  (0.01,  0.1,  1,  and  2)  using 
the  Nicolis  [1992]  model  with  stochastic  forcing  q 2  =  0.2.  The  contour  plot  covers  the  half  domain  due  to  £0  >  £noise. 


with  the  initial  error  £0  bounded  by 

Snoise  '  SO  i  “■ 


dj\ 

dto 


=  0, 


di2 

dto 


0,  £o 


noise  • 


(27) 


Furthermore,  equations  (18)  and  (19)  become  ordinary  differential 
equations 


It  d2TX  _ 

«o-Sio)^o+  —  ^T- 


«o  -  gi  o) 


d'^2  .  q2i 0  <P"l2 


dto 


+ 


^0 


=  -2ti 


(24) 


(25) 


4.3.  Analytical  Solutions 

[18]  Analytical  solutions  of  equations  (24)  and  (25)  with  the 
boundary  conditions  (26)  and  (27)  are 

l 

—  —  2  f  — 2zg 

Tl(£,0,  Snoise,  e)  =  /  T  ?eXP(— T) 

r  J  <r 

U 


with  the  boundary  conditions, 


y 

/ 


“exp  {—^-x)dx 

r 


dy 


(28) 


Tj  =  0, 


T2  =  0,  £0  =  s; 


(26) 
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Figure  3.  Contour  plots  of  t2  (£,0 ,  E,„oise  >  e)  versus  $,0 ,  E,noise)  for  four  different  values  of  £  (0.0 1 , 0. 1 ,  1 ,  and  2)  using  the 
Nicolis  [1992]  model  with  stochastic  forcing  q 2  =  0.2.  The  contour  plot  covers  the  half  domain  due  to  (;0  >  £,noise. 


l 


Unoise 


where  Q,  =  We  and  ^noise  =  ^.noise/6  are  the  nondimensional 
initial  condition  error  and  noise  level  scaled  by  the  tolerance  level 
e,  respectively.  For  given  tolerance  and  noise  levels  (or  user  input) 
the  mean  and  variance  of  VPP  can  be  calculated  using 
equations  (28)  and  (29). 

4.4.  Dependence  of  iq  and  t2  on  (£„,  lnoise/<z) 

[19]  To  investigate  the  sensitivity  of  t2  and  t2  to  Q^noise 
and  If,  the  same  values  are  used  for  the  parameters  in  the 


stochastic  dynamical  system  (22)  as  were  used  by  Nicolis 
[1992]: 

a  =  0.64,  g  =  0.3,  q-  =  0.2.  (30) 

Figures  2  and  3  show  the  contour  plots  of  ti (^0,^noise,  e)  and 
T2(^0>  incise,  e)  versus  (£,o,£,noise)  for  f°ur  different  values  of  e 
(0.01,  0.1,  1,  and  2).  The  following  features  can  be  obtained:  (1) 
For  given  values  of  (QuCnoise)  the  same  location  in  the 

contour  plots)  both  n  and  t2  increase  with  the  tolerance  level  e. 
(2)  For  a  given  value  of  tolerance  level  e  both  iq  and  t2  are 
almost  independent  on  the  noise  level  ^noise  (contours  are  almost 
parallel  to  the  horizontal  axis)  when  the  initial  error  Q,  is  much 
larger  than  the  noise  level  £,noiSe-  This  indicates  that  the  effect  of 
the  noise  level  S,noise  on  iq  and  t2  becomes  evident  only  when 
the  initial  error  S;0  is  close  to  the  noise  level  £,nojse.  (3)  For  given 
values  of  (e,  S;noise)  both  tj  and  t2  decrease  with  increasing 
initial  error  Q). 
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Figure  4.  Dependence  of  ti  (S;0,  E;noise,  e)  on  the  initial  condition  error  Q,  for  four  different  values  of  £  (0.01,  0.1,  1, 
and  2)  and  four  different  values  of  random  noise  £,no;se  (0.1,  0.2,  0.4,  and  0.6)  using  the  Nicolis  [1992]  model  with 
stochastic  forcing  q 2  =  0.2. 


[20]  Figures  4  and  5  show  the  curve  plots  of  D  (QnCnoisei e) 
and  T2(F,0,£,noise5£)  versus  £,0  for  four  different  values  of  toler¬ 
ance  level  s  (0.01,  0.1,  1,  and  2)  and  four  different  values  of 
random  noise  S;noise  (0.1,  0.2,  0.4,  and  0.6).  The  following 
features  are  obtained:  (1)  t,  and  t2  decrease  with  increasing 
Q),  which  implies  that  the  higher  the  initial  error,  the  lower  the 
predictability  (or  VPP);  (2)  T!  and  t2  decrease  with  increasing 
noise  level  f)noise,  which  implies  that  the  higher  the  noise  level, 
the  lower  the  predictability  (or  VPP);  and  (3)  ti  and  t2  increase 
with  the  increasing  e,  which  implies  that  the  higher  the  tolerance 
level,  the  longer  the  VPP.  Note  that  the  results  presented  in  this 
section  are  for  a  given  value  of  stochastic  forcing  ( q 2  =  0.2) 
only. 

4.5.  Dependence  of  Ti  and  t2  on  Stochastic  Forcing  q2 

[21]  To  investigate  the  sensitivity  of  Ti  and  t2  to  the  strength 
of  the  stochastic  forcing  q2,  we  use  the  same  values  for  the 


parameters  (a  =  0.64  and  g  =  0.3)  in  equation  (30)  as  were  used 
by  Nicolis  [1992],  except  q2,  which  takes  values  of  0.1,  0.2,  and 
0.4^Figures  6  and  7  show  the  curve  plots  of  Ti(Q},Cnoise)92)  and 
'U^oiS.noise!1?2)  versus  £,0  for  two  tolerance  levels  e  (O.land  1), 
two  noise  levels  S;noise  (0.1  and  0.6),  and  three  different  values  of 
q2( 0.1,  0.2,  and  0.4)  representing  weak,  normal,  and  strong 
stochastic  forcing.  Two  regimes  are  found:  (1)  ti  and  t2  decrease 
with  increasing  q2  for  large  noise  levels  (^noise  =  0.6),  (b)  x,  and 
t2  increase  with  increasing  q2  for  small  noise  levels 
(Xnoise  =  0.1),  and  (3)  both  relationships  (increase  and  decrease 
of  Ti  and  t2  with  increasing  q2  are  independent  of  s.  These 
indicate  the  existence  of  stabilizing  and  destabilizing  regimes  of 
the  dynamical  system  depending  on  stochastic  forcing.  For  a 
small  noise  level  the  stochastic  forcing  stabilizes  the  dynamical 
system  and  increases  the  mean  VPP.  For  a  large  noise  level  the 
stochastic  forcing  destabilizes  the  dynamical  system  and 
decreases  the  mean  VPP. 
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Figure  5.  Dependence  of  ^noise> e)  on  the  initial  condition  error  Q,  for  four  different  values  of  e  (0.01,  0.1,  1, 
and  2)  and  four  different  values  of  random  noise  £,no;se  (0.1,  0.2,  0.4,  and  0.6)  using  the  Nicolis  [1992]  model  with 
stochastic  forcing  q2  =  0.2. 


[22]  The  two  regimes  can  be  identified  analytically  for  small  is  negligible.  The  solution  (31)  becomes 
tolerance  level  (e  — ►  0).  The  initial  error  Q,  should  also  be  small 
(Q)  e).  The  solution  (28)  becomes 


JimTi  (Q),  Snoise’  2) 


Nnoisc’  -) 

E — >0 


(33) 


1  tr 
2a  -  q2 


1 


(31) 


The  Lyapunov  exponent  is  identified  as  a  —  q2/ 2  for  dynamical 
system  (22)  [Has  ’minskii,  1980],  Fora  small  noise  level  (E;noise  =  1) 
the  second  term  in  the  bracket  of  the  right-hand  of  equation  (31), 


R  =  - 


V  Fr1 

2a  -  q2  ^nolse 


(32) 


which  shows  that  the  stochastic  forcing  (q  ^  0)  reduces  the 
Lyapunov  exponent  (ct  —  q2/ 2),  stabilizes  the  dynamical  system 
(22),  and  in  turn,  increases  the  mean  VPP.  On  the  other  hand,  the 
initial  error  Q,  reduces  the  mean  VPP. 

[23]  For  a  large  noise  level  £(no;se  the  second  tern  in  the  bracket 
of  the  right-hand  side  of  equation  (31)  is  not  negligible.  For  a 
positive  Lyapunov  exponent,  2a  —  q2  >  0,  this  term  is  always 
negative  (see  equation  (32)).  The  absolute  value  of  R  increases 
with  increasing  q2  (remember  that  S;noise  <  1  and  Q,  <  1).  Thus  the 
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Figure  6.  Dependence  of  ti  (S,0,  £,noisei  ?2)  on  the  initial  condition  error  Q,  for  three  different  values  of  the  stochastic 
forcing  q2  (0. 1 ,  0.2,  and  0.4)  using  Nicolis  model  with  two  different  values  of  e  (0. 1  and  1 )  and  two  different  values  of 
noise  level  ^noise  (0.1  and  0.6). 


term  R  destabilizes  the  one-dimensional  stochastic  dynamical 
system  (22)  and  reduces  the  mean  VPP. 


5.  Conclusions 

1.  The  valid  prediction  period  (a  single  scalar)  represents  the 
model  predictability  skill.  It  depends  not  only  on  the  instantaneous 
error  growth  but  also  on  the  noise  level,  the  tolerance  level,  and  the 
initial  error.  A  theoretical  framework  was  developed  in  this  study 
to  determine  the  mean  (iq)  and  variability  (t2)  of  the  valid 
prediction  period  for  a  nonlinear  stochastic  dynamical  system.  The 
probability  density  function  of  the  valid  prediction  period  satisfies 
the  backward  Fokker-Planck  equation.  After  this  equation  it  is  easy 
to  obtain  the  ensemble  mean  and  variance  of  the  valid  prediction 
period. 

2.  Uncertainty  in  ocean  (or  atmospheric)  models  is  caused  by 
measurement  errors  (initial  and/or  boundary  condition  errors), 
model  discretization,  and  uncertain  model  parameters.  This 


provides  motivation  for  including  stochastic  forcing  in  ocean 
(atmospheric)  models.  The  backward  Fokker-Planck  equation  can 
be  used  for  evaluation  of  ocean  (or  atmospheric)  model 
predictability  through  calculating  the  mean  valid  prediction 
period. 

3.  For  an  autonomous  dynamical  system,  time-independent 
second-order  linear  differential  equations  are  derived  for  iq  and  t2 
with  given  boundary  conditions.  This  is  a  well-posed  problem,  and 
the  solutions  are  easily  obtained. 

4.  For  the  Nicolis  [1992]  model  the  second-order  ordinary 
differential  equations  of  iq  and  t2  have  analytical  solutions, 
which  clearly  show  the  following  features:  (1)  decrease  of  iq  and 
t2  with  increasing  initial  condition  error  (or  with  increasing 
random  noise)  and  (2)  increase  of  iq  and  t2  with  increasing 
tolerance  level  e. 

5.  Both  stabilizing  and  destabilizing  regimes  are  found  in  the 
Nicolis  [1992]  model  depending  on  the  stochastic  forcing.  For  a 
small  noise  level  the  stochastic  forcing  stabilizes  the  dynamical 
system  and  increases  the  mean  VPP.  For  a  large  noise  level  the 
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Figure  7.  Dependence  of  T2(£o,  ^noisei  q2)  on  the  initial  condition  error  S;0  for  three  different  values  of  the  stochastic 
forcing  q 2  (0.1,  0.2,  and  0.4)  using  the  Nicolis  [1992]  model  with  two  different  values  of  e  (0.1  and  1)  and  two 
different  values  of  noise  level  !;nois(,  (0.1  and  0.6). 


stochastic  forcing  destabilizes  the  dynamical  system  and  decreases 
the  mean  VPP. 
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